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According to Lovelock’s theorem, the Hilbert-Einstein and the Lovelock actions are indistinguishable 
from their field equations. However, they have different scalar-tensor counterparts, which correspond to the 
Brans-Dicke and the Lovelock-Brans-Dicke (LBD) gravities, respectively. In this paper the LBD model of 
alternative gravity with the Lagrangian density .^lbd = li; [4> (^ + is developed, 

where "RR and respectively denote the topological Chern-Pontryagin and Gauss-Bonnet invariants. The 
field equation, the kinematical and dynamical wave equations, and the constraint from energy-momentum 
conservation are all derived. It is shown that, the LBD gravity reduces to general relativity in the limit tUL ^ 
unless the “topological balance condition” holds, and in vacuum it can be conformally transformed into 
the dynamical Chern-Simons gravity and the generalized Gauss-Bonnet dark energy with Horndeski-like or 
Galileon-like kinetics. Moreover, the LBD gravity allows for the late-time cosmic acceleration without dark 
energy. Finally, the LBD gravity is generalized into the Lovelock-scalar-tensor gravity, and its equivalence 
to fourth-order modified gravities is established. It is also emphasized that the standard expressions for the 
contributions of generalized Gauss-Bonnet dependence can be further simplified. 
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I. INTRODUCTION 


As an alternative to the various models of dark energy with 
large negative pressure that violates the standard energy con¬ 
ditions, the accelerated expansion of the Universe has inspired 
the reconsideration of relativistic gravity and modifications of 
general relativity (GR), which can explain the cosmic accel¬ 
eration and reconstruct the entire expansion history without 
dark energy. 

Such alternative and modified gravities actually encode the 
possible ways to go beyond Lovelock’s theorem and its nec¬ 
essary conditions [1], which limit the second-order field equa¬ 
tion in four dimensions to H- Ag^y = 

i.e. Einstein’s equation supplemented by the cosmological 
constant A. These directions can allow for, for example, 
fourth and even higher order gravitational field equations [2- 
5], more than four spacetime dimensions [6, 7], extensions of 
pure pseudo-Riemannian geometry and metric gravity [7, 8], 
extra physical degrees of freedom [9-12], and nonminimal 
curvature-matter couplings [13, 14]. From a variational ap¬ 
proach, these violations manifest themselves as different mod¬ 
ifications of the Hilbert-Einstein action, such as extra curva¬ 
ture invariants, scalar fields, and non-Riemannian geometric 
variables. 

For the Lovelock action in Lovelock’s theorem and the 
Hilbert-Einstein-A action, it is well known that they yield the 
same field equation and thus are indistinguishable by their 
gravitational effects. When reconsidering Lovelock’s theo¬ 
rem, we cannot help but ask whether the effects of these two 
actions are really the same in all possible aspects. Is there any 
way for the two topological sources in the Lovelock action to 
show nontrivial consequences? As a possible answer to this 


question, we propose the Lovelock-Brans-Dicke gravity. 

This paper is organized as follows. In Sec. II, the Lovelock- 
Brans-Dicke gravity is introduced based on Lovelock’s the¬ 
orem, and its gravitational and wave equations are derived 
in Sec. III. Section IV studies the behaviors at the infinite- 
Lovelock-parameter limit wl —> oa, and Sec. V derives the 
constraint from energy-momentum conservation. Section VI 
shows that in vacuum the Lovelock-Brans-Dicke gravity can 
be conformally transformed into the dynamical Chern-Simon 
gravity and the generalized Gauss-Bonnet dark energy with 
Horndeski-like or Galileon-like kinetics. Then the possibility 
of realizing the acceleration phase for the late-time Universe 
is discussed in Sec. VII. Finally, in Sec. VIII the Lovelock- 
Brans-Dicke theory is extended to the Lovelock-scalar-tensor 
gravity, and its equivalence to fourth-order modified gravities 
is analyzed. Throughout this paper, we adopt the sign conven¬ 
tions = dyYl^ - ■ and R^y = 

with the metric signature (-, -H ■¥ -h). 


II. LOVELOCK-BRANS-DICKE ACTION 


An algebraic Riemannian invariant R — R [gap , Rafipv) in 

the action J d^x s/^R generally leads to fourth-order gravi¬ 
tational field equations by the variational derivative 

Sgr- Sggy "afAg"’) ” ' 


Lovelock found out that in four dimensions the most general 
action leading to second-order field equations is [1] 
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where A is the cosmological constant, {a , b] are dimensional 
coupling constants, and without any loss of generality we have 
set the coefficient of R equal to one. Also, refers to the 
totally antisymmetric Levi-Civita pseudotensor with eoi 23 = 

obtained from 

each other by raising or lowering the indices with the metric 
tensor. In Eq.(2), and Q respectively refer to 

the Chern-Pontryagin density and the Gauss-Bonnet invariant, 
with 


g R^- AR^,R>^^ + R^aviiR'^"''''. ( 3 ) 

The variational derivatives dieapfivR^'’ !Sg^'' and 

yield total derivatives which serve as boundary 
terms in varying the full action Eq.(2). The Chern-Pontryagin 
scalar eap/jyR^'' is proportional to the divergence of the 

topological Chern-Simons four-current [11], 

with 

and similarly, the topological current for the Gauss-Bonnet in¬ 
variant is (see Refs. [15, 16] for earlier discussion and Ref. [17] 
for further clarification) 


= -dX with 


Hence, the covariant densities eap^yR^'' in 

Eq.(2) make no contribution to the field equation dSjdg^'' - 0. 

According to Lovelock’s theorem, one cannot tell whether 
Einstein’s equation R^y - jRg^iv = comes from the 

customary Hilbert-Einstein action 

‘She = d^xyPgR + S,n, (6) 


or from the induced Lovelock action 


‘Sl = J" +‘S„ 


with 
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[11, 18, 19]. In Eqs.(2), (6) and (7), the matter action Sm 
is given in terms of the matter Lagrangian density by 
S„, - Jd^x and the stress-energy-momentum den¬ 

sity tensor is defined in the usual way by [20] 

= d‘^xyPgTj^^6g^'' with 

-2 

■ V=5 ■ 


The indistinguishability between >Sl and .She from their 
field equations begs the question: Does Einstein’s equation 
come from .Sl or .She? E there any way to discriminate them? 

Recall that GR from .She has a fundamental scalar-tensor 
counterpart, the Brans-Dicke gravity [9], 


‘Sbd - 


J" 
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BD 




+ Sm with 


-2bd - -r^ ^R -— Va([)V' 
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which proves to be a successful alternative to GR that passes 
all typical GR tests [21], and it is related to GR by 


-SfuE - 
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That is to say, Brans-Dicke firstly replaces the matter-gravity 
coupling constant G with a pointwise scalar field ([)(x“) in ac¬ 
cordance with the spirit of Mach’s principle, G 4)“\ and 
further adds to the action a formally canonical kinetic term 
VQ,(j)V“([) governing the kinetics of tt)(x“). Applying 

4* 

this prescription to the Lovelock action Eq.(7), we obtain 


-2lbd = 
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where the Lovelock parameter wl E a dimensionless constant. 
Based on Eq.(12), we obtain what we dub as the Lovelock- 
Bmns-Dicke (henceforth LBD) gravity with the action 


where for simplicity we switch to the denotation 

^RR (8) 

for the Chern-Pontryagin density, as the symbol *RR has 
been widely used in the literature of Chern-Simons gravity 


' Note that not to confuse with the more common “Lovelock action” for 
the topological generalizations of the Hilbert-Einstein action to generic N 
dimensions that still preserves second-order field equations, as in Ref. [6]. 


‘^LBD - J" d'^X yXg -S^LBD + Sm , (13) 

or the Lanczos-Lovelock-Brans-Dicke gravity, as Lovelock’s 
theorem is based on Lanczos’ discovery that an isolated *RR 
or g in the action does not affect the field equation [22]. 

Unlike the 6(*RR)/6gf^'' and 6(^Xgg)|6g^^'' in 6Svl6gi^\ the 
[(\>dCRR)] I6g>^'' and [q^b(V=?^)] |6g^^'' for b.SLBo/bg^*' are no 
longer pure divergences, because the scalar field c])(x“) as a 
nontrivial coefficient will be absorbed into the variations of 
*RR and yj-^g when integrating by parts. Hence, although 
.Sl and .She are indistinguishable, their respective scalar- 
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tensor counterparts >Slbd and >Sbd are different. 

Note that the cosmological-constant term -2A in Eq.(2) is 
temporarily abandoned in otherwise, it would add an ex¬ 
tra term -2Ac|) to ^lbd, which serves as a simplest linear po¬ 
tential. This is primarily for a better analogy between the LED 
and the Brans-Dicke gravities, as the latter in its standard form 
does not contain a potential term y(4)), and an unspecified po¬ 
tential y(4)) would cause too much arbitrariness to ^lbd- 

Also, Lovelock’s original action Eq.(2) concentrates on the 
algebraic curvature invariants; in fact, one can further add 
to Eq.(2) the relevant differential terms nf?, n*RR, and a0 
(□ = denoting the covariant d’Alembertian), while 

the field equation will remain unchanged. This way, the grav¬ 
itational Lagrangian density in Eq.(7) is enriched into 

1 / q*RR \ 

- \R +——+ bg + caR + da*RR + eag], (14) 

16;rG \ ) 


III.l. Coupling to the Chern-Pontryagin invariant 

The Chern-Pontryagin density *RR in J^bbd measures the 
gravitational effects of parity violation through J d^x<i^*RR 
for its dependence on the Levi-Civita pseudotensor. In ad¬ 
dition to Eq.(8), *RR is related to the left dual of the Riemann 
tensor via 

*RR = \ {eap,vR^\s)R‘'^^' = *Ral3y6R°^^^' ■ (17) 

Applying the Ricci decomposition = Capy6 + \{gayRi}6- 

Sa^Rpy Ep^Ray ~ 8PyRa&^ ~ 6 {^EayEpS ~ EadSPy'^ R Eq.(17) 
and using the cyclic identity Capys + Caysp + Caspy - 0 for the 
traceless Weyl tensor, one could find the equivalence 

*RR = *CC := i {eap^.yC^\s) = *CapyeC‘^'^^^ , (18) 


with {c, d, e] being constants, and its Brans-Dicke-type coun¬ 
terpart extends Eq.(12) into 


.if = 


1671 


R + ——*RR + bg\- — V„(|3V“(1) 
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-I- (|5 + dD*RR + eag ] 


(15) 


where cf) ■ (cnR + dD*RR + eag) have nontrivial contribu¬ 
tions to the field equation. However, unlike *RR and yFsg 
which are divergences of their respective topological cur¬ 
rent as in Eqs.(4) and (5), {dR, □*RR, ag] are total deriva¬ 
tives simply because the d’Alembertian □ satisfies yf^aQ - 
da ( y[—EE“^dp^ when acting on an arbitrary scalar field 0; 
in this sense, these differential boundary terms which contain 
fourth-order derivatives of the metric are less interesting than 
*RR and g. In this paper, we will focus on the LBD grav¬ 
ity .ifuBD Eq.(12) built upon the original Lovelock action and 
Lovelock’s theorem, rather than Eq.(15) out of the modified 
action Eq.(14). 


III. GRAVITATIONAL AND WAVE EQUATIONS 


In this section we will work out the gravitational field equa¬ 
tion ^.Slbd/^.?^’" = 0 and the wave equation (5 .Slbd/<^<)> = 0 
for the LBD gravity. Lirst of all, with 5gap - -EafiEpy^E^'’^ 
^^ip ^ {'^adEa-p + ^pdEa-a “ ^tr^Eap), and the Palatini 
identity = '^p (^ryo.) - [23], for the first term 

(pR in .jSfLBD it is easy to work out that 

1 (5(V=^(1)R) 1 , . 

(16) 

where = means equality by neglecting all total-derivative 
terms which are boundary terms for the action. 


which indicates that the Chern-Pontryagin density is confor¬ 
mally invariant [15] under a rescaling Q.{xF ' Ei^y of 

the metric tensor. 

With the Chern-Simons topological current in Eq.(4), 
one can integrate by parts and obtain J d^x (j)*RR - 

-4 f d'^x ([) (5^/:^) = -4 / d'^xd^ + 4 fd^x (5^c[.) R/^. 
Hence, instead of directly varying (t)*RR with respect to the 
inverse metric, we firstly vary the four-current R^ by the Levi- 
Civita connection. It follows that 

6 J" d'^x cf)*RR - J" ('^7'4*) dR'^ 

= ijd^x (d^(i>) (V^dgay - Vydga^) 

= - 2 J /x [{d^<^) 

^-4jd^x [(5^(|r) d“^PyVpR; + dg,y (19) 

- 4 ^d\ [(5^c^) e,apyFR\ + (20) 

where, in the third row we expanded (JL^^ and made use of 
the cancelation F''^^Vadg^y = 0 due to the skew-symmetry 
for the indices ^ v; in the fourth row, we applied the re¬ 
placement pRy - '^jRp in accordance with the 

relation 

'^“Raiipy - ^pR/jy “ yRjip, ( 21 ) 

which is an implication of the second Bianchi identity ^jRafipy 
+'^yRaiiyp + ^pRafiyy - 0; in the last step, we raised the indices 
of dgay to 6g°’'' and thus had the overall minus sign dropped. 
In Eq.(20) we adopted the usual notation cl^c[) s and 

note that *R^/y * R^opy since in general the 

metric tensor does not commute with partial derivatives and 
thus d^d^<!s^ - g^^dfi ^ g^^g^^dfid^c)). Relabel the in- 
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dices of Eq.(19) and we obtain the variational derivative 


1 6ii>*RR) 
^F8 ^g^'' 


and 


( 22 ) 

+ Idadu^ ■ [*R\^ + *R"/^) . 

Compared with Eq.(16), does not contain a -^(^*RRgi^y 
term, because *RR by itself already serves as a covariant den¬ 
sity as opposed to the usual form for other curvature 

invariants. 

Note that the nonminimal coupling between a scalar field 
and *RR is crucial to the Chern-Simons gravity; however, 
its original proposal Ref. [11] had adopted the opposite ge¬ 
ometric system which uses the metric signature (-H,-), 

the conventions {R “''' ’ Ein¬ 

stein’s equation - ^Rg/jv — -87TGTj^\ and the defini¬ 
tion *RR = -*RapysR'^^^^ = -\{eap^,vR^\s)R‘'^^^- This has 
caused quite a few mistakes in the subsequent Chern-Simons 
literature that adopt different conventions, and we hope the de¬ 
tails in this subsection could correct these misunderstandings. 
Also, in Eq. (22), the quantities K^, *RR, R^y, *Rpfiav} 

have the same values in both sets of sign conventions. See our 
note Ref. [24] for further clarification of this issue. 


III.2. Coupling to the Gauss-Bonnet invariant 


The third term c])^ in J^lbd represents the nonminimal cou¬ 
pling between the scalar field and the Gauss-Bonnet invariant 
0 - R^ - AR^ -I- where we have employed the straight¬ 
forward abbreviations := Rq./?R“^ and R^, := Ra/jpvR"^^'' 
to denote the Ricci and Riemann tensor squares. Eollowing 
the standard procedures of variational derivative as before in 
we have 


5 (V^4>^) 
~v/=f^ 




SgMV 


Sgt^y 


Sgi^y 


(23) 


with 


5(^R2) 

5g^^y 

5(^R2) 

5g^y 

<5(q^R2) 

Sg^y 


= 24)RR^v + 2 (g^vD - V^Vy) {(^R) (24) 

= 2<\,R^ “R„v + □ {<\>R^y) - V„V, (c^R/) 

-V„V^ (^R/) + g^yV^Vp (^R"^) (25) 

= 2^R^„^^R, + 4V^V“ (q^R„^^v) , (26) 


where total-derivative terms have been removed. Recall that 
besides Eq.(21), the second Bianchi identity also has the fol¬ 
lowing implications which transform the derivative of a high- 


rank curvature tensor into that of lower-rank tensors plus non¬ 
linear algebraic terms; 

^^Rap = \ ^pR (27) 

V^V“R„^ = 1 dR (28) 

^^^‘'Ra^pv = aR^v - ^V^V,R Ra^pyR^^ - R^'^R^v (29) 

V“V^R^ + V"V,R,^ = V^V,R - 2R,,pyR‘^P + 2R/R^(30) 

Using Eq.(21) and Eqs.(27)-(30) to expand the second-order 
covariant derivatives in Eqs.(24)-(26), and putting them back 
into Eq.(23), we obtain 


1 5{^F-g^9) 
yFg Sg^'' 


with 


(2RR^, - 4R/R„, - 4R„^^,R“^ + 2R^„^^R/^’') 
+27^ — V^Vyjc]) H-4 R^"Vq,Vv(]) h- 4Ry“VQ,V^([) 

-4R^yn([i - 4g^yR“^VQ.V^(|) -h 4RQ.^^yV^V“c|) 

-^(i>9giiy, (31) 

where the second-order derivatives (□, VqVv, etc] only act on 
the scalar field ([). 

However, we realize that Eq.(31) is still not the ul¬ 
timate expression. In four dimensions, is pro¬ 

portional to the Euler-Poincare topological density, Q - 
{{e.pyiR-^^'^^) ■ and the in¬ 

tegral J dx^ y[^9 equates the Euler characteristic x{M) 
of the spacetime. Thus f dx^ y/^9 = 32n'^-^%{M) = 
0. Based on Eqs.(24)-(26), one could easily obtain the 
Bach-Lanczos identity from the explicit variational derivative 
s{^9)l6g^‘\ 

2RR^y-4R^^R,y-AR,^pyR">^ + 2R,,pyRy‘^^^ = ^gg^y, (32) 

with which Eq.(31) can be best simplified into 

77;^®) =2R (g^yD - V^Vy) + 4R/V„Vy^ + 4Ry “V„V^ 

-4R^yD(i> - 4g^y ■ R^^V^V^c^ + 4R,^^yV^V>, (33) 
whose trace is 

gPV7^(GB) ^ 2Rn^ - 4R“^V„V^^. (34) 

In the existent literature, the effects of the generalized and 
thus nontrivial Gauss-Bonnet dependence for the field equa¬ 
tions are generally depicted in the form analogous to Eq.(31), 
such as the string-inspired Gauss-Bonnet effective dark en¬ 
ergy [12] with ^ = j^R-^df.ipd^'if-Viip) +f{tp)g, as well 
as the R -H f{g) [3], the f(R,g) [4] and the f(R,g,Xn) [14] 
generalized Gauss-Bonnet gravities. Here we emphasize that 
the Gauss-Bonnet effects therein could all be simplified into 
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the form of Eq.(33). 


IV. THE Wl ^ oo LIMIT AND GR 


III.3. Gravitational field equation 


Collecting the results in Eqs.(16), (22), and (33), we finally 
obtain the gravitational field equation 

cl) iR^,v - ^ ^V^c|5Vv(|t - 

+ - V^V,) 

(CP) 

where vanishes for all spherically symmetric or confor¬ 
mal flat spacetimes. Eq.(35) yields the trace equation 

-^R+ — V„<l)V“ct) H- (3 -I- 2bR) net) - 4M"^V„V«cl) = 8kT^’^\ 

(36) 

where is always traceless, = 0 - this is not a 

surprise because it equivalently traces back to the effects of 
the dual square *CC of the traceless Weyl tensor. 

Note that in existent studies the invariants *RR and Q have 
demonstrated their importance in various aspects. Eor exam¬ 
ple, as shown by Eq.(6) of Ref. [25] [recall the equivalence 
*RR - *CC in Eq.(18)], in the effective field theory for the 
initial cosmic inflation, the only leading-order fluctuations to 
the standard inflation action in the tensor modes are the parity- 
violation Chern-Pontryagin and the topological Gauss-Bonnet 
effects. 


III.4. Wave equations 


Straightforward extremization of <Slbd with respect to the 
scalar field yields the kmematical wave equation 

= + (37) 

with □([) = We regard Eq. (37) as “kine- 

matical” because it does not explicitly relate the propagation 
of c[) to the matter distribution .if„, or 

Combine Eq.(37) with the gravitational trace equation (36), 
and it follows that 

(2wl + 3-1- 2bR) □([)=- {-^’‘RR + bg\ t|) 

\ yFs I (38) 

-H 87rE'”> H- 

which serves as the generalized Klein-Gordon equation that 
governs the dynamics of the scalar field. 


Erom the dynamical equation (38), we obtain 

□d) =-i-1 -i-^'‘RR + bg\<h 

^ 2wl + 3+2M\ 

+ 4M‘’^V„V/5(t)|. (39) 

The topology-gravity coupling strengths {a,b} should take 
finite values - just like the Newtonian constant G for 
matter-gravity coupling. Similarly the curvature invariants 
{R, *RR, g) for a physical spacetime should be finite, and we 
further assume the scalar field ([) to be nonsingular. Thus, in 
the limit ojl Eq.(39) yields □([) = O and 

(|) = (([)) +©(—) = ^-I-©(—), (40) 

\ wl / G \ wl / 

where (([)) denotes the expectation value of the scalar field and 
we expect it to be the inverse of the Newtonian constant 1/G. 
Under the behaviors Eq.(40) in the infinite wl limit, we have 
= 0 = and the field equation (35) reduces to 

become Einstein’s equation - jRgfiv - 

On the other hand, from Eq.(39) we can also observe that 
□([) s 0 in the special situation 


- 4fcR“^V<,V^([. -H ^ *RR + bg'^(^ = , (41) 

and the scalar field becomes undeterminable from the dynam¬ 
ical equation (39). 

The term -4bR°'^Vjj(^ comes from the trace g^''H^^^\ 
while *RR and g are respectively related to the topological 
instanton number [15] and the Euler characteristic. Thus, all 
terms on the left hand side of Eq.(41) are related to topologi¬ 
cal effects nonminimally coupled with cj), and they cancel out 
the trace of the matter tensor. In this sense, we call Eq.(41) 
the topological balance condition. 

Putting □<[) s 0 and the condition Eq.(41) back into the 
trace equation (36), we obtain 

^V„(|)V> = R + -^*RR + bg (42) 

4) ^Fg 

- -wlD In (]), (43) 

where in the second step we further made use of the expansion 
□ Incj) = V“(iV„(|t) = -iV„cl)V“(|) inc]) = 
for □([) s 0. Thus it follows that 


wlV„ (In ([)) V“ (In ^ R + -^*RR + bg . (44) 

Fs 

Eor cul ^ oo, this equation gives the estimate 


||V„(ln(» 


R + -^*RR + bg 


OJl 


O 



(45) 
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which integrates to yield Incf) = constant + Hence, 

(j) satisfies 


(() ~ tjto + (46) 

where the constant 4>o is the average value of tj). In ac¬ 
cordance with Eq.(42) and the estimate Eq.(46), the term 
(v^4)Vy(l) - in the field equation (35), 

which arises from the source -^VQ,4)V“(t) in .Sled, will not 
vanish. This way, the wl ^ oo limit could not recover Ein¬ 
stein’s equation and GR in situations where the topological 
balance condition Eq.(41) holds, although the existence of 
such solutions remains to be carefully checked. 

This is similar to the Brans-Dicke theory given by the action 
Eq.(lO), which recovers GR in the limit wbd ^ unless the 
stress-energy-momentum tensor has a vanishing trace - 
0 [26], such as the matter content being radiation with Rrad = 
|pi-ad and = -Prad + 3Prad = 0. 

V. ENERGY-MOMENTUM CONSERVATION 


Moreover, for the scalar field, we have 


V^' 


■ Jr + - 


(50) 


where the kinematical wave equation (37) has been employed. 

Eor the Chern-Pontryagin and the Gauss-Bonnet parts in 
Eq.(35), consider the componential actions .Scp = f d‘^x^*RR 
and >Sgb = f d‘^x Under an arbitrary infinitesimal 

coordinate transformation xf^ ^ xf^ + 6x^, where 6x^ - is 
an infinitesimal vector field which vanishes on the boundary, 
so that the spacetime manifold is mapped onto itself. Then 
>Scp and .Sgb vary by 


SScp = -J d"^x(pdf,(^^*RR) = J d''x*RR(d,j(t>)^^, (51) 
dScB = - J d'^xclyd^ J d^x . 


(52) 


In modified gravities with the generic Lagrangian density 
^ = f{R, %, ■■■), where % = Riiga/i , Raiifiv , '^yRaii/Sv , • • • , 
VyjVy, ... '^y„RaiJi 3 v) and the “■ ■ ■ ” in = / refer to arbi¬ 
trary curvature invariants beyond the Ricci scalar, the energy- 
momentum conservation is naturally guaranteed by Noether’s 
law or the generalized contracted Bianchi identities [27] 


1 <5[V=g/(R,R„---)] ^ 


which can be expanded into 


(47) 


(48) 


where /« df(R, %, ■ ■ ■ )/5R, f% := df(R, ■ ■ ■ )/&Ri, and 
R-Ijv - {fKid%) 15g^''. However, in the more generic situa¬ 
tions of scalar-tensor-type gravities with ^ - /(<]), R, %, • • ■) 
-i-iD-(c[), VQ.4)V“t])) where nonminimal couplings between the 
scalar fields and the curvature invariants are involved, such as 
the LBD proposal under discussion, the conservation problem 
is more complicated than pure tensorial gravity. 

Now let’s get back to the LBD field equation (35). By the 
coordinate invariance or the diffeomorphism invariance of the 
matter action Sm in which is neither coupled with the 
curvature invariants nor the scalar field cf), naturally we have 
the energy-momentum conservation = 0 for the mat¬ 

ter content. Thus, the covariant derivative of the left hand 
side of Eq.(35) should also vanish. With the Bianchi identity 

[R/jv - jRgiJv) - 0 and the third-order-derivative commu¬ 
tator (VyD - nVy) (]) = -R^yV^c]), it follows that 










■ (.g^vD - V^Vy) C|) 


= --RVy(^. (49) 


Eor the first step in Eqs.(5 1) and Eqs.(52), one should note that 
H- is a particle/active transformation, under which 
the dynamical tensor fields transform, while the background 
scalar field tt)(x“) and the coordinate system parameterizing 
the spacetime manifold remain unchanged [28]. On the other 
hand, the inverse metric transforms by g'^'' + 6g^'' with 

6 g^^'' = = V^^'' -I- and thus we have 

dScP^ 

= -2 (53) 

SSgb^ 2 

= -2 (54) 

Comparing Eqs.(51) with (53), and Eqs.(52) with (54), we 
obtain the relations 


= (55) 

Adding up Eqs.(49), (50), (55), and (56), one could find 
that the covariance divergence for the left hand side of the 
field equation (35) vanishes, which confirms the energy- 
momentum conservation in the LBD gravity. 

(CP^ 

Eqs.(55) and (56) for the nontrivial divergences of H^jy 
and H'^^\ by their derivation process, reflect the breakdown 
of diffeomorphism invariance for >Scp and >Sgb in 'Slbd- They 
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have clearly shown the influences of nonminimal cjj-topology 
couplings to the covariant conservation, as opposed to the 
straightforward generalized Bianchi identities 




/ 1 6*RR\ 

\ ^Fs Sg^'' I 


= 0 and 


lJ_shFg^\ 

I ^Fg ^8^'' I 


= 0. 


(57) 


VI. CONFORMAL TRANSFORMATIONS 


The standard LBD action >Slbd in Eq.(12) can be trans¬ 
formed into different representations by conformal rescaling 
of the spacetime line element, which geometrically preserves 
the angles between spacetime vectors and physically retains 
local causality structures. 


VI. I. Dynamical Chern-Simons gravity 


As a simplest example, consider the specialized <Slbd in 
vacuum and for spacetimes of negligible gravitational effects 
from the nonminimally (|)-coupled Gauss-Bonnet term. With 
S,„ = 0 and b -0, Eq.(12) reduces to become 




16;r 


f 


d^x 




Wl, 


V„(|)V“(|5 -H a^*RR 


(58) 


Eor a pointwise scaling held Q = G(x") > 0, we can rescale 
the metric of the original frame into g^y via 
it follows that gf,y = g^'' = ^F8 = 

and2 [10] 

R^n^[R + 6n(ln Q) - 6|“^5„(ln Q) «9^(ln Q)]. (59) 

Hence, for the reduced LBD action Eq.(58), the conformal 
transformation 


1 


(60) 


along with the redefinition of the scalar held {?7 = &(x"), (|) = 
c[)(i?)) lead to 


.5 = 


IbTrG 


/ 


ctx 




R- 


2(Ul + 3 
2cl)(l?)2 





+a (^idYRR 


( 61 ) 


where the scalar held ■& no longer directly couples to the 
Ricci scalar R, and thus the 6Ci(lnQ) component in Eq.(59) 
has been removed as it simply yields a boundary term 
6 f da [y^(9“(ln Q)j d^x for the action. Also, Eq.(61) has uti¬ 


lized the fact that the (l,3)-type Weyl tensor and thus 

*RR - *CC - *CC - *RR are conformally invariant. It is 
straightforward to observe from Eq.(61) that the kinetics of i? 
is canonical for wl > -3/2, noncanonical for wl < -3/2, and 
nondynamical for ujl = -3/2; here we are interested in the 
canonical case. Eor the specialization 

dd = ± V2wl + 3 —, (62) 

which integrates to yield 

= +V2wl + 3 In-^, (63) 

([30 


where ([to is an integration constant, or inversely 


c[) = 4)0 exp + 


V2a)L + 3 / 
the action Eq.(65) Anally becomes 




- f 

16;rG J 


dF ^(R-iv^dV^d 


+a4>o exp + 




V2u)l + 3 


*RR 


(64) 


(65) 


Hence, the conformal rescaling g^y - f^y/G4) along with the 
new scalar held i7(x“) recast the reduced LBD action Eq.(58) 
into Eq.(65), which is an action for the dynamical Chern- 
Simons gravity [19], though the nonminimal d-*RR coupling 
is slightly more complicated than the straightforward d*RR as 
in the popular Chern-Simons literature. Moreover, the confor¬ 
mal invariance of *RR guarantees that the effect of J d^x (p*RR 
could never be removed by conformal transformations. 

Note that the matter action SmigiivFm) would be trans¬ 
formed into S,n{giAvlG^,\l/,„) (in general S does not contain 
derivatives of the metric tensor [20]), which are different in 
the 4>-‘5m or 4>--2]« couplings; consequently fails to be 
conformally invariant unless it is traceless = 0 [10]. This 
is why we focus on the vacuum situation. 


VI.2. Generalized Gauss-Bonnet dark energy 


Similarly, in vacuum and for spacetimes of negligible 
Chern-Simons parity-violation effect, <Slbd reduces into 








( 66 ) 


Under the local rescaling gfjy i-» g^^y - Gl?g^y for the metric, 
the Gauss-Bonnet scalar satisfies [29] 


^ Compared with - £1^ [s + 6Qn/f2 llg'^^daCidpO.ICi?-^, Eq.(59) best coetBcient of R is reset into unity. Moreover, by employing InH instead of 

isolates pure-divergence terms and thus most simplifies the action once the the transformations R^R becomes skew-symmetric \oR^R. 
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g = - 8^“^V<,V;3(ln Q) - 8^“^V<,(ln Q)V;3(ln Q) + AR aOn Q) + 8 [nan Q)]^ - 8a(ln O) ■ V„(ln Q)V"(ln Q) 

- 8V„V/5(lnO) ■ V“V^(lnQ) - 16V„V;3(lnQ) ■ V“(lnQ)V^(lnQ)l. 


(67) 


Set the factor of conformal transformation to be Q = -^Gcf) 
so that the Ricci scalar decouples from the scalar field, and 
redefine the scalar field via cl)(x“) i,c>(x“) = V2wl + 3 In ^ 

or equivalently (jj = cjDo exp f j; then it follows that 


In Q = - In ( 
2 

1 


+ -InG 

2 


[f) 1 1 

-—== + -Incfio + - InG. 
2 '^2 (Ul + 3 2 2 


( 68 ) 


With {In 4)0, In G) being constants, substitution of Eq.(68) into 
Eq.(67) yields 




7C(V^) = - 2R‘^^ {iVjlpip + 

+ Hip ■ {iR + 2nip — (70) 

- 2Vj!0 ■ (V“V^^ + V“(^V^(,o) . 

Here one can observe that since the coefficient Q in - 
Q exactly neutralizes the O'* in Eq.(67), the nonmin- 

imally 4>-coupled Gauss-Bonnet effect J d'^x could 

never be canceled by a conformal rescaling g/jy - 
Hence, the reduced LBD action Eq.(66) is finally transformed 
into 


/ \ 1 

which generalizes the canonical Gauss-Bonnet dark energy in 
vacuum S = ^ f {r - -F f{ip)g) [12] 

by the Horndeski-like [30] or Galileon-like [31] kinetics in 
'K{V(p) for the scalar field. 

Note that in the two examples just above, because of the 
nonminimal coupling to the scalar field 4)(x“), negligible 
Gauss-Bonnet effect does not imply a zero Euler character¬ 
istic x(^) = 3 ^/ - 0 for the spacetime, and 

similarly, negligibility of the Chern-Simons effect does not 
indicate a vanishing instanton number J *RR d'^x - 0, either. 

Also, for the actions of the Chern-Simons gravity and the 
Gauss-Bonnet dark energy in the Jordan frame, in which a 
scalar field is respectively coupled to *RR and g, we cannot 


help but ask that why the scalar field is not simultaneously 
coupled to the Ricci scalar? We have previously seen from 
Eq.(47) that all algebraic and differential Riemannian invari¬ 
ants stand equal in front of the generalized Bianchi identities, 
so are there any good reasons for the scalar field to discrim¬ 
inate among different curvature invariants? We hope that the 
LBD gravity help release this tension (at least in empty space- 
times), as the scalar field cf) indiscriminately couples to all the 
LBD invariants {R,*RR,g], and the LBD gravity takes the 
Chern-Simons gravity and the Gauss-Bonnet dark energy as 
its reduced representations in the Einstein frame. 


VII. COSMOLOGICAL APPLICATIONS 


Having extensively discussed the theoretical structures of 
the LBD gravity, in this section we will apply this theory to the 
Eriedman-Robertson-Walker (ERW) Universe and investigate 
the possibility to realize the late-time cosmic acceleration. 


VII.l. Generalized Friedmann and Klein-Gordon equations 

The field equation (35) can be recast into a GR form, 

R,v - \Rg,. = Fg,s , (72) 

where F - 8;7r, and Ggff =1/4) denotes the effective gravita¬ 
tional coupling strength. -H F^y + =: 

compromises the total effective stress-energy-momentum ten¬ 
sor, with 

^and (73) 

/c-r* = (V;,Vv - 4> + ^ |v^4>Vv4> - ^.g/<vVff4>v“4)j. 

Note that besides the effects of the source term -^V(j,4>V“4> 
in .5fLBD via 6 (- ^Vff4>V“'4') the (v^Vy - g^yo) 4) 

part from I6g^'' is also packed into . More¬ 

over, with the four distinct components of sharing the 
same gravitational strength 1/4), Eq.(72) implicitly respects 
the equivalence principle that the gravitational interaction is 
independent of the internal structures and compositions of a 
test body or self-gravitating object [21]. 

Eor the ERW metric of the flat Universe with a vanishing 
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spatial curvature index. 


and for pressureless dust = 0, it integrates to yield 


ds^ - —dt^ + a(t)^ ^ 

i=\ 




{m) 

^0 . 


-3p 


( 81 ) 


*RR - 0 due to the maximal spatial symmetry, while the Ricci 
and Gauss-Bonnet scalars are respectively 


ad + iP' / . 

Q = 24^ = 2AH^ (h + H^) , 


(75) 


where we have assumed a power-law scale factor 

a - aoP with p > 1. (82) 

Here {ao, P) are constants, and P > 1 so that d > 0. Similarly, 
we also take a power-law ansatz for the scalar field, 

^ (83) 


where overdot denotes the derivative over the cosmic comov¬ 
ing time, and H dja represents the time-dependent Hubble 
parameter. Thus, an accelerated/decelerated flat Universe has 
a positive/negative Euler-Poincare topological density. With 
a perfect-fluid form = diag [-p, P, P, P] assumed for each 
component in [in consistency with the metric signature 
(-,-1- -I- -I-)], the cosmic expansion satisfies the generalized 
Friedmann equations 


Based on Eqs.(81)-(83), the dynamical wave equation (38) 
with = -p,„ for dust yields 


.2 „(in) 


Y ( 2 wl + 3 ) ( 3 P - 1 + y) = 

cpoao 

p 3 (P- 1 ) , ( 3 P -3 -hy) 

+24b'—^^—- - 


(84) 




(76) 


a 

a 


1 f 2c(j 

— \k^ (p„, + 3P„) + 3$ + 

6(|3 [ ([. 

-H \2bH^'^ + 12b[2H -I- H^) I, 


(77) 


where - 0 for FRW. Moreover, the kinematical wave 
equation (37) and the dynamical wave equation (38) respec¬ 
tively lead to 


2a>\ /.. . \ ad + cP 


wl 


• ([) H- 24bH^ 


(78) 


and in the late-time (large f) approximation it reduces to 


1 2 (m) 

y(y + 3 P- 1 )^^ 


([>0 (2a>L + 3) 


which can be satisfied by 


Y = 2 - 3p and c()o = 


2 (m) 


fl3(2wL + 3)(2-3p)' 


Moreover, the first Friedmann equation (76) leads to 


J’”) 




12b 


P'y 


(85) 


( 86 ) 


(87) 


a^cljo ■ ' 2 ■ f- 

and with Eq.(86), in the late-time approximation it becomes 


/ ad + d^\/- ■\ id 

\2col + 3 -h 12b j(cf) H- 3H^) = 24bH^-<\> 

( d ■■ ad + 2cP ■ \ 

-([) H- -T-//([) . 

a a^ I 


(79) 


3p2 = (2 - 3P) (2 wl + 3) - 3P(2 - 3P) + ^(2 - 3P)^ (88) 

For p = 2, Eq.(88) trivially holds for an arbitrary wl, while 
for p 4 2, we have P in terms of wl via 


In principle, one could understand the evolutions of the scale 
factor a{t) and the homogeneous scalar held cj)(f) by (prob¬ 
ably numerically) solving Eqs.(76)-(79). However the solu¬ 
tions will be complicated, so we will start with some solution 
ansatz for {a{t), c[)(f)), which are easier to work with. 


VII.2. Cosmic acceleration in the late-time approximation 


2(fa;L+ 1) 
3 wl + 4 


(89) 


Note that Eqs.(86) and (89) reuire wl 4 -4/3, wl 4 -3/2 
(p 4 2), and /? 4= 2/3; they are simply consequences of 
the power-law-solution ansatz and the late-time approxima¬ 
tions rather than universal constraints on a>L, and according 
to Eq.(89), the last condition [} + 2/3 trivially holds with 
P ^ 213 for wl ^ As a consistency test, the kinemat¬ 
ical equation (37) yields 


The physical matter satisfies the continuity equation 

p,„ + 3H{p„, + P„,) = 0, (80) 


WLI - r + 3ySrI = 3^2^ - i) + i2b 




(90) 
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with the late-time approximation 


and attributes to the effect of 


COL Qr - r + 3/3r) = 3/3(2/? - l), (91) 

which holds for Eqs.(86) and (89). Substituting Eqs.(81), 
(82), (83), (86) and (89) into the second Eriedmann equation 
(77), we obtain 


a _ 2(tOL + 1 )(col + 2) ^_2 
a (3wl + 4)2 

and the deceleration parameter reads 

— o-Peff\ u;l + 2 

2\ peffj 2 (wl + 1)' 


(92) 


(93) 


=2R - V^V,)/3(^) - AR^yOM,^) + 

4/?/V„Vv/3((^) + 4/?/V„V^/3((^) - (97) 

It is straightforward to derive the kinematical wave equation 
by (5.SLST/b4) = 0, which along with the trace of Eq.(95) could 
yield the dynamical wave equation, and they generalize the 
wave equations (37, 38) in the LED gravity. The wave equa¬ 
tions however will not be listed here as the interest of this 
section is only the field equation SSjSg^'" - 0. 


Eqs.(92) and (93) clearly indicate that the late-time acceler¬ 
ation could be realized for-2 < wl < -1 (col + -4/3, VIII.2. Equivalence of LST with fourth-order gravities 

wl + -3/2), although this domain of col makes the kinetics 

of the scalar field noncanonical. . „ , . . ^ n . . . . 

It IS well known that the j{R) gravity is equivalent to the 

nondynamical (i.e. wbd = 0) Brans-Dicke gravity [27], and 

such equivalence holds for the LED gravity as well. Consider 

the fourth-order modified gravity 


VIII. LOVELOCK-SCALAR-TENSOR GRAVITY 

.if = 


VIII.l. Erom LED to Lovelock-scalar-tensor gravity 

for which the field equation is 


16;rG 


f(R,g) + h 


I ’•RR \ 


(98) 


The LBD gravity can be generalized into the Lovelock- 
scalar-tensor (LST) gravity with the action 

‘5lst = J" yPg^LST + S,n and 

■^LST = l^j^Q (94) 

V(c»), 

where {//((])), cu(([))} are generic functions of the scalar field, 
and y(4)) is the self-interaction potential. Note that this time 
Newton’s constant G is included in the overall coefficient 
l/16;7rG of -SfLsx. as is the case of the ordinary scalar-tensor 
gravity. The gravitational field equation is 

/i(4>) - v^v^)/i(([)) 

- ^ (95) 

+ ^^V(<i,)g,y + 

where W^y denotes the contribution from f 2 {<^yRR, 

=25^/2(^) ■ [e^,ay^‘^R^y + 

+ 2d,dpf2{<^) ■ {*R%^ + 


IrRiiv + (g/jvD - - 2^f{R,Q)gfiv 

+ <//(CP) + .^(GB) ^ 8^p(m) 

^ fiv * ^ /iV ’ 

where 

+ 2d,diih.RR ■ [*R%^ + , 

and 

^2R [g,yU - V,Vy)fg - AR^yOfg + 

43?/V„V,/p + 4/?/V„V^/p- 

with /r = fR(R,g) = df(R,g)ldR, fg = fgiR.Q) = 
df{R,@)ld@, snA h.RR = dhCRRI 2 f^)ldCRRI Eor 

the nondynamical LST gravity with 01(4)) = 0 in Eq.(95), 
compare it with Eq.(99) and at the level of the gravitational 
equation, one could find the equivalence 

/Kct*) = /r > /3(4') = fg > /2(4') = h-RR, 2') 

v{i,)^-mg) + fRR. 

In the y(cf)) relation we have applied the replacement /i (cf)) = 
/«, and note that y(4)) does not contain a fgff term which 
has been removed from because of the Bach-Lanczos 

identity Eq.(32). 


(99) 


( 100 ) 


( 101 ) 
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VIII.3. Partial equivalence for “multi-scalar LED gravity” 


Removing the oil term in Eq.(35) and then comparing it 
with Eq.(99), one could find that an equivalence between the 
nondynamical LED gravity (now equipped with an extra po¬ 
tential -t/(cl)) in ^lbd) and the f{R,Q) + gravity 

would require = cj) = Ii.rr, and {/(c])) = -f{R,@) 

+fRR. These conditions are so restrictive that the f{R,0) + 

gravity would totally lose its generality. Instead, in¬ 
troduce three auxiliary fields {xi,X2,X3) and consider the dy¬ 
namically equivalent action 


+ /^(T3) + 


fiX\^X2) + fx, ■iR-X\)+fx 2 -iS-Xi) 
*RR 


2 F 8 


-Xi 


"t" 5 


(103) 


its variation with respect to x\, Xi^ and X 2 separately yields 
the constraints 


^^,(R-^i) = o, fx^,(g-x2)-o, 

^ , (*RR ) n (104) 

and -Xi) ^ 0, 

where/^, df{xuX 2 )ldxj, fxjxj - d^fiX\^X 2 )ldx], 

dh(x3)ldxi := d^f(x2)ldxl. If fxixi^ fxixi ^^^>^ 2 x 2 

do not vanish identically, Eq.(104) leads to;^fl - R, X 2 - G 
and ^3 = Redefining the fields {xi’X2,X3} by 


<>^fx,’ M' = fx2’ ( 105 ) 

and setting 

* R.R 

E(cl), \\f, cp) =c|) ■ R(^, v) -F V ■ ^(cj), v) H- cp ■ —=(cp) 

yFg 

( *RR \ 

-^(cp)|, (106) 

then the f{R, @) + h j gravity is partially equivalent to the 

following “multi-scalar LED gravity” carrying three nondy¬ 
namical scalar fields 

-S? = Y^(ct)^ + cp-^+¥^-(107) 

where the coupling coefficients {a, b} appearing in .ifuBo have 
been absorbed into the scalar fields {cp, \|/}. Also, by “partially 
equivalent” we mean that Eq.(106) as is stands is only par¬ 
tially on-shell; to recover Eq.(98) from the multi-field action 
ofEq.(107), one would have to add extra Lagrange multipliers 
identifying the different fields, but this would break the exact 
equivalence between such modified Eq.(107) and Eq.(98). 


IX. CONCLUSIONS AND DISCUSSION 


The Hilbert-Einstein action .She ancl the Lovelock action 
.Sl yield identical field equations and thus are observationally 
indistinguishable. However, the former takes the Erans-Dicke 
gravity as its scalar-tensor counterpart, while the latter’s com¬ 
panion is the LED gravity, and these two theories are different. 

We have extensively studied the theoretical structures of 
the LED gravity, including the gravitational and wave equa¬ 
tions, the ordinary wl —> °o limit that recovers GR, the un¬ 
usual wl ^ 00 limit satisfying the topology balance condition 
Eq.(41) and thus departing from GR, the energy-momentum 
conservation, the conformal transformations into the dynam¬ 
ical Chern-Simons gravity and the generalized Gauss-Eonnet 
dark energy, as well as the extensions to LST gravity with its 
equivalence to fourth-order modified gravity. 

We have taken the opportunity of deriving the field equa¬ 
tion to look deeper into the properties of the Chern-Pontryagin 
and Gauss-Eonnet topological invariants. Especially, for the 
Gauss-Eonnet dark energy as well as the f{R, G) and 
f{R,G,-^m) gravities, the contributions of the generalized 
Gauss-Eonnet dependence could be simplified from the pop¬ 
ular form like Eq.(31) into our form like Eq.(33). 

An important goal of alternative and modified gravities is 
to explain the accelerated expansion of the Universe, and we 
have applied the LED theory to this problem, too. It turned 
out that the acceleration could be realized for -2 < cul < -1 
under our solution ansatz. Note that our estimate of cosmic 
acceleration in Sec. VII is not satisfactory. Lor example, the 
kinematical equation (90) clearly shows that because of the 
higher-order time derivative terms arising from the ^G depen¬ 
dence, the simplest solution ansatz {(jt = 4)of^, a - with 
{P=constant, y^constant) are not compatible with each other 
unless the late-time approximation is imposed, while such ap¬ 
proximations further lead to the behaviors analogous to the 
Erans-Dicke cosmology [32]. 

Section VII has shown that, the effects from the parity- 
violating Chern-Pontryagin term ^*RR are ineffective for 
the LRW cosmology because of its spatial homogeneity and 
isotropy. However, it is believed that ^*RR could have 
detectable consequences on leptogenesis and gravitational 
waves in the initial inflation epoch [33] where c]) acts as the 
inflaton field. The inflation problem usually works with the 
slow-roll approximations c]) « (jt <s // and requires the exis¬ 
tence of a potential well V(4)); thus, at least for the description 
of the initial inflation, the LED gravity should be generalized 
to carry a potential: 


.^LBD - 


1 

16;r 


cl) 


^R + Cl 


*RR 



^v„(^v>-y(^) , 

(108) 


with y(4)) = 2Ac|) being the simplest possibility. 

Our prospective studies aim to construct the complete his¬ 
tory of cosmic expansion in LED gravity [probably equipped 
with y(4))], throughout the dominance of radiation, dust, and 
effective dark energy. Moreover, it is well known that primor¬ 
dial gravitational waves can trace back to the Planck era of 
the Universe and serve as one of the most practical and effi- 
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cient tests for modified gravities, so it is very useful to find 
out whether the gravitational-wave polarizations carry differ¬ 
ent intensities in this gravity. There are also some other prob¬ 
lems from the LBD gravity attracting our attention, such as its 
relation to the low-energy effective string theory. We will look 
for the answers in future. 
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